Relaxation electrodynamics of superinsulators

Superinsulators offer a unique laboratory realizing strong interaction phenomena like confinement and asymptotic freedom in quantum materials. Recent experiments evidenced that superinsulators are the mirror-twins of superconductors with reversed electric and magnetic field effects. Cooper pairs and Cooper holes in the superinsulator are confined into neutral electric pions by electric strings, with the Cooper pairs playing the role of quarks. Here we report the non-equilibrium relaxation of the electric pions in superinsulating films. We find that the time delay \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t_{\mathrm {sh}}$$\end{document}tsh of the current passage in the superinsulator is related to the applied voltage V via the power law, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t_{\mathrm {sh}}\propto (V-V_{\mathrm p})^{-\mu }$$\end{document}tsh∝(V-Vp)-μ, where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V_{\mathrm p}$$\end{document}Vp is the effective threshold voltage. Two distinct critical exponents, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu =1/2$$\end{document}μ=1/2 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu =3/4$$\end{document}μ=3/4, correspond to jumps from the electric Meissner state to the mixed state and to the superinsulating resistive state with broken charge confinement, respectively. The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu =1/2$$\end{document}μ=1/2 value establishes a direct experimental evidence for the electric strings’ linear potential confining the charges of opposite signs in the electric Meissner state and effectively rules out disorder-induced localization as a mechanism for superinsulation. We further report the memory effects and their corresponding dynamic critical exponents arising upon the sudden reversal of the applied voltage. Our observations open routes for exploring fundamental strong interaction charge confinement via desktop experiments.

www.nature.com/scientificreports/ The two critical voltages appear as two kinks in the I(V) characteristics 1 . However in the low-resistive samples near the SIT, V c1 and V c2 may often appear very near to each other, and experiments performed in this region can fail to separate them. The corresponding observation of only one kink has led to a speculation 24,25 that the superinsulating behavior might arise due to disorder-induced-localization 26 . This idea originated from the partial explanation of the shape of the jump in the I(V) curves by charge overheating 24 , which was erroneously taken as a model of the nature of the superinsulating state, while it has only to do with one possible mechanism by which this state is dynamically destroyed. Here we report the first ever experiment on the dynamical behaviour of superinsulators and we show that, even when the two critical voltages are so close that they cannot be distinguished in the I(V) curves, they are clearly identified by their different dynamic critical exponents relating the delay in the current response to a sudden voltage pulse, t sh ∝ (V − V p ) −µ , where V p is the effective threshold voltage. Furthermore, we show that the exponent µ = 1/2 describing the jump from the Meissner state to the mixed state is a direct confirmation of the linear binding potential between charges caused by electric strings in the confined state, while the second exponent µ = 3/4 reflects the quantum phase slips in the 1D edge channel of the intermediate bosonic topological insulator phase separating superinsulation from superconductivity. These results fully support the monopole-induced confinement nature of superinsulators and strongly disfavour the assumption that disorder-induced-localization may model the superinsulating state.

Experiment
The measurements are taken on the superinsulating NbTiN films grown by the atomic layer deposition (ALD) technique. The low-temperature experiments are carried out in a 3 He/ 4 He dilution refrigerator Triton 400 equipped with the measuring RF lines. To measure the current when the low frequency, 2 × 10 4 Hz, pulsed voltage is applied, the measuring circuit consisting of an Agilent 81150A pulse generator, a Femto HCA highspeed current amplifier with a bandwidth of 0 -2 MHz, and a Tektronix DPO 7354C oscilloscope, is used. The time of the voltage rise/fall in this circuit does not exceed 500 ns. Figure 1a displays a typical response of the current in a superinsulating film to a square pulse voltage, as compared to the dynamic response of the normal insulating state at higher temperatures. Contrary to the normal insulator, in which the current increases continuously to the steady state, in the superinsulator at the first moment, no measurable current value appears in the system. Then, after some delay time, t sh , a measurable current arises and quickly reaches saturation. When the voltage is turned off, the current decreases exponentially. These features of the current response are observed provided the pulse duration is at least 1.5 times longer than the time t sh . In Figure 1b we present the dependencies of the shift time of the rising edge t sh upon the amplitude V of the pulse signal for several samples which differ from each other in resistivity and hence in the resulting threshold voltages. Note that with the increase in resistivity, the increase in the threshold voltage is observed. As the threshold voltage grows, the rising edge shift time first increases, then reaches a maximum, and then, finally, begins to decrease. We find that the shift time t sh is a power law function of the reduced voltage V − V p , i.e., t sh ∝ (V − V p ) −µ , where V p is a fitting parameter which depends only upon the sample. At some voltage, the power µ abruptly changes. In addition, at voltages less than V c2 , the effect of shifting the front edge of the pulse and measurable current level are observed only when the signal is periodic with the period being less than 1 ms. The small AC component promotes the depinning of the charges at the ends of electric flux tubes in an exact www.nature.com/scientificreports/ analogy to the AC current component promoting the vortex mobility and the resistive behaviour in the mixed state of superconductors in the experiments on superconductors 23 . The voltage V c2 corresponds to the second threshold voltage of the film 1 , see Fig. 2. At the same time, the voltage V p is smaller than the threshold voltage, V c2 , which has a hysteretic behavior presented in Fig. 2. As we will show in the next section, V p must be identified with the first threshold voltage V c1 . Figure 3 shows similar results for voltage jumps from +V to −V . Unlike the above 0 → |V | jumps, now as the sign of the voltage changes, the current also instantly changes its sign. However, in the further time interval, t pm , the current decreases and, having reached its minimum depending on the voltage value, again rapidly increases to saturation. Thus, when the sign of the voltage changes, the "initial" dynamic state is first destroyed and only then a new state is formed. The time t pm is also a power function of the reduced voltage V − V p , with the same probe voltage but with another exponent. Note that, when the sign of the voltage changes, there are no features associated with the threshold voltage V c2 . At the same time, at voltages V between V p and V c2 the requirement for a signal to be periodic with a period of less than 1 ms is persisting.   www.nature.com/scientificreports/ screened Coulomb interactions and bound by the linear potential V (x 1 , x 2 ) = σ |x 1 − x 2 | , where σ is the string tension of the confining electric string and x i are the coordinates of the two particles. The critical value of the potential when current starts to pass is V c1 = σ L , with L being the sample length, corresponding to a configuration in which the two opposite charges reach the sample boundaries.

Results
In center of mass coordinates we can reduce the problem to one particle of the reduced mass m/2 at the origin, subject to the potential energy U(r) = 2eσ r , with r denoting the separation distance. The binding energy of the ± charge pair is 2eV c1 ≈ 0.3 eV, where we used the average value of the threshold voltages in Fig. 1. This is far smaller than the reduced mass 0.5 MeV of a Cooper pair. We can thus safely use classical considerations to compute the behaviour of the centre of mass coordinate. A positive linear potential results in an attractive constant force F a = 2eσ = 2eV c1 /L . Let us now apply an external potential V. This amounts to an additional repulsive force F r = 2eV /L . The total acceleration of separation is thus a = (2/m)F tot = (4e/mL)(V − V c1 ) . The corresponding equation of motion separating the charge-hole system is Current starts to pass when the charges of opposite sign reach the sample boundaries, i.e. when r(t cr ) = L . The delay for current passage is thus given by inverting Eq. (1), This shows that the observed value µ = 1/2 of the dynamical critical exponent for the current passage is a direct consequence of the linear binding potential between charges of the opposite sign and that V p must be identified with the first critical voltage V c1 , thereby fully confirming the confinement nature of superinsulators. Note, that even if the I(V) measurements are not accurate enough to disentangle the two critical voltages, they can still be clearly identified by their two different dynamic critical exponents.

Quantum dissipation.
For higher applied voltages the system jumps directly over the second critical voltage V c2 beyond which the superinsulation is completely destroyed. Had the resulting state at temperature T = 20 mK, been an Ohmic resistive state, where it is the velocity of charges which is proportional to the applied force rather than the charges' acceleration, the above simple calculation shows that we would have expected the critical exponent µ = 1 , which is clearly not the case. Instead, for NbTiN, this state is a BTI [27][28][29] , with a topological gap for bulk excitations and charge transport mediated by the symmetry-protected edge modes 18,30 . Only at higher temperatures one would have observed a thermally activated bulk conductance. At the lowest temperatures, we expect that the quantum dissipation occurs only due to quantum phase slips 31-33 across the 1D edge conducting channels.
The temperature-independent but velocity-dependent friction is, indeed, a fundamental characteristic of the 1D bosonic conduction channels in the quantum phase slips regime 34 . Such a friction has recently been experimentally confirmed, although there is not yet full quantitative agreement with theory 35 . The velocity dependent damping coefficient implies an anomalous relation between the applied force F ∝ (V − V cr ) and the velocity, with K being the Luttinger parameter of the bosons, see 36 for a review. This is the same power law that has been originally obtained in 37 , where it was shown that the quantum phase slip nucleation rate Ŵ plays the role of the applied force in the quantum wire dissipation. Repeating the above calculation we obtain an arrival time scaling which matches the experimental result for K = 5/3.

Switching delay.
Let us now address the delay t ± observed upon the sudden switching the direction of the applied electric field. Suppose that we start off in a situation with the negative voltage having the magnitude larger than V c2 and therefore, the corresponding steady-state current goes in the negative direction. When switching to the positive voltage of the same magnitude, the system will momentarily pass through the mixed and Meissner states of the superinsulator. In this process, the strings form, which are initially taut in the negative direction. But when the applied field passes through the zero, the charges can flow freely in the opposite direction for a while, until the string completely "turns around", i.e., until the tension builds up in the opposite direction. This explains the initial bump in the current, which then decreases when the string tension increases to oppose the inverted applied force. Since the applied voltage is larger than V c2 , at a certain point the mixed state is again destroyed to form a steady current, now flowing in the positive direction. The delay t pm is thus another, albeit qualitative, confirmation of the string-confinement picture of the superinsulation. A quantitative derivation of the observed scaling of t pm with the power µ = 2/3 is beyond the present theoretical framework.

Discussion and conclusion
The results of our measurements of dynamical exponents describing the relaxation of the current passage confirm the existence of two electric states of superinsulators, the Meissner state, where the electric field is expelled from the sample and the mixed state where electric field penetrates superinsulator in form of the electric filaments (2) t cr = www.nature.com/scientificreports/ or electric strings. These two states of a superinsulator constitute an exact duality to the conventional Meissner and mixed states in a superconductor. Furthermore, the exponent µ = 1/2 for the transition from the Meissner to the mixed state provides a direct validation of the linear binding potential created by the electric strings between charges of opposite sign in the electric Meissner state. Our results rule out any superinsulation model not containing this linear binding potential preventing charge transport at voltages not exceeding the critical voltage. In particular, it strongly disfavors disorder-induced localization as a mechanism for current suppression and hyperactivation.

Data availability
The data that support the findings of this study are available from the corresponding authors upon reasonable request.